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ABSTRACT 

Aims. A proper determination of the abundance gradient in the Milky Way requires the observation of objects at large galactiocentric 
distances. With this aim, we are exploring the planetary nebula population towards the Galactic Anticentre. In this article, the discovery 
and physico-chemical study of a new planetary nebula towards the Anticentre direction, IPHASXJ05 253 1.1 9+28 1945.1 (PNG 178.1- 
04.0), is presented. 

Methods. The planetary nebula was discovered from the IPHAS survey. Long-slit follow-up spectroscopy was carried out to confirm 
its planetary nebula nature and to calculate its physical and chemical characteristics. 

Results. The newly discovered planetary nebula turned out to be located at a very large galactocentric distance (Doc = 20.8±3.8 kpc), 
larger than any previously known planetary nebula with measured abundances. Its relatively high oxygen abundance (12+log(0/H) = 
8.36±0.03) supports a flattening of the Galactic abundance gradient at large galactocentric distances rather than a linearly decreasing 
gradient. 
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1. Introduction 

The INT Photo metric H-Alpha Survey of the Northern Galactic 
Plane (IPHAS. iDrew et alj|2005 iGonzalez-Solares et al.ll2008h 
is significantly increa sing the number of det ected Galactic plan- 
etary nebulae (PNe). IViironen et al.l d2009a) published a list of 
781 new IPHAS PN candidates, whereas the discovery and 
detailed study of fiv e new IPHAS PNe were presented by 
IViironen et ai1(l2009bl) . 

PNe towards the Galactic Anticentre are of special interest as 
this is where the objects with the largest galactocentric distances 
are most easily found and studied. Finding new distant PNe will 
make a crucial contribution to the knowledge of the Galactic 
abundance gradient, since their rich emission line spectra can 
yield accurate abundances. The abundance gradient is a topic 
that is still under debate (see e.g. lMaciel & Costall2009l and ref- 
erences therein). In addition to the shortage of known objects at 
large galactocentric distances, the main problem is th e difficulty 
of reliably measuring PN distances (see Kwok 2000, for exam- 
ple). The new faint (and therefore potentially distant) IPHAS 
Anticentre PN candi dates combined with new methods fo r dis- 
tance determination (ISale et al.ll2009t; iFrew & Parkerll2006b can 
make a significant contribution to the problem. 

The location of the newly discovered PN IPHASX 
J052531. 19+281945.1, at only four degrees from the direction 
of the Galactic Anticentre, motivated us to study this object fur- 
ther. In Sect. [2] the photometric and spectroscopic observations 
are presented; in Sect. [3]a physical and chemical analysis is car- 
ried out; in Sect. |4]its distance is studied; in Sect. [5]a discussion 
is provided, and in Sect.[6]the main conclusions are presented. 

2. Observations and data reduction 

The PN was discovered by visually inspecting the IPHAS Ha - 
r' mosaics of the Anticen tre region. The semi-au tomated search- 
ing method described in IViironen et al.l (1200 9 a) missed the ob- 
ject due to its relatively large size (~ 10")- The IAU approved 
name for this IPHAS object is, after its coordinates, IPHASX 
J052531. 19+281945.1 (X for extended) whereas the IAU com- 
mon PN name would be PNG 178.1-04.0. To abbreviate, in the 
following we will use the acronym IACPN (IPHAS Anti-Centre 
Planetary Nebula). 

2.1. IPHAS imaging 

The IPHAS observations of IACPN were carried out on 
November 2003 using the 2.5m Isaac Newton Telescope and its 
Wide Field Camera under seeing condition of 1.3" FWHM in the 
Ha images. Fig.[T]presents the Ha, r 1 and i' images of the object. 
The exposure times are 120 s for Ha and 10 s for r' and i' im- 
ages. More detail s about IPHAS obs erva tions and the data reduc- 
tion a re given in IDrew et ail (120051) and Gonzalez-Solares et alJ 
(2008). 

IACPN shows a roughly spherical morphology. It is hardly 
visible in the IPHAS i' filter (see Fig. [T}. No objects are listed in 
SIMBAD within almost 5' from the position of IACPN. 

2.2. SPM 2.1 -m spectroscopy 

The first spectroscopic observations of IACPN were carried out 
at the Obs ervatorio de San Pedro M artir (Mexico), and were pre- 
sented by Mampaso et al. (2005) in a conference paper. Four 
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20 min exposures in the red and five 20 min exposures in the 
blue were obtained through a narrow (1") slit using the 2.1-m 
telescope and the B&Ch spectrograph (2 A pix -1 dispersion, 3 
A resolution, from 3720 to 7470 A). In addition, two 10 min 
exposures in the blue were obtained with a very wide slit in or- 
der to include the whole nebula and measure its total Hp flux. 
Furthermore, the MEZCAL echelle spectrograph was used with 
the same telescope to obtain a high dispersion spectrum (0.1 
A pix 1 ) around the Ha- line. 

The MEZCAL spectrum allowed to measure the expansion 
and radial velocities of the nebula: the ~ 10" diameter shell has 
a radial velocity of v^r — 13.5 km s -1 and expands at 17 
km s -1 in the [Nn] 6583 A line (M ampaso et al.ll2005l) . From 
the B&Ch spectra the PN nature of the nebula was confirmed 
and its absolute H/3 flux was measured, giving F(H/3) = 3.34 x 
10~ 14 erg cm" 2 s -1 . Preliminary chemical abundances were pre- 
sented by Mampaso et al] d2005l) . 

2.3. WHT + ISIS spectroscopy 

The William Herschel Telescope (WHT) + ISIS spectrograph 
observations of IACPN were carried out in service mode on 
April 6, 2007. The seeing was ~ 0.9" and the weather condi- 
tions were good, except for the presence of light cirrus during 
the second half of the night. Three exposures of 20 min of the 
target were obtained in both the blue (3570-51 15 A; dispersion 
0.8 A pix" 1 , resolution 3.2 A) and the red (5540-10500 A, 1.8 
A pix -1 dispersion, 6.5 A resolution) arms with the slit located 
at PA. = 79°, i.e. at the parallactic angle corresponding to the 
moment of observation. The slit width was 1". In addition, bias, 
lamp flats, and arc exposures were obtained. Unfortunately, no 
suitable observations of spectrophotometric standards are avail- 
able for that night, as all standard stars were observed with a 
different instrumental configuration than our target. For t his rea - 
son, a spectrum of the standard star BD+332642 dOkd [l990). 
obtained on April 2, 2007 with the same setup as for the nebula, 
was used to flux calibrate its spectrum. Nevertheless, the blue 
and red observations of this standard star were obtained at dif- 
ferent times of the night, and therefore under different weather 
conditions. This prevents a precise matching of the blue and red 
sides of the WHT spectrum, and therefore, an accurate determi- 
nation of Ay using the HafH/3 flux ratio. For this reason, addi- 
tional spectra of the nebula and standard star were acquired as 
described below. Finally, and given that there is calibration data 
available for BD+33 2642 only up to 9200 A dOkdll990l) . the 
fluxes for lines at longer wavelengths will not be considered. 
The data reduction and calibration were carried out using the 
standard IRAF routines for longslit spectroscopy. 

2.4. NOT + ALFOSC spectroscopy 

A good measurement of the reddening of the nebula, Ay, is im- 
portant for its distance determination (see Sec. |4) and additional 
spectra of IACPN were secured in service mode on September 
6, 2010 using the ALFOSC instrument at the Nordic Optical 
Telescope. These spectra were obtained so as to cover the wave- 
length range from ~3800 A to 6800 A, i.e. to include simulta- 
neously the main Balmer lines (H6, Hy, H/3, and Ha). The night 
was photometric and the seeing was 0.7". Two 20 min expo- 
sures were obtained. The spectral dispersion was 1.5 A pix 1 
and resolution 6 A. The slit was positioned according to the par- 
allactic angle at the moment of observation (PA. = 102°) and 
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Fig. 1. IPHAS Ha, r' and i' discovery images of IACPN, from left to right. The image sizes are 1' x 1', North is up and East to the 
left. 
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Fig. 2. Spatial profiles extracted from the NOT spectrum for rep- 
resentative lines from ions with different ionisation potentials. 

the slit width used was 1". Bias, lamp flats, and arc exposures 
were obtained, and the standard star G191-B2B was observed. 
The data reduction and calibration were carried out using the 
standard IRAF routines for longslit spectroscopy. 



3. Spectroscopic analysis 

The spatial distribution of the emission from several lines in the 
NOT spectrum is presented in Fig. [2] and shows that IACPN 
exhibits a strong ionisation stratification, with higher excitation 
emissions coming from regions progressively closer to the cen- 
tre. For this reason, care has to be exercised when comparing 
line fluxes from the WHT and NOT spectra. The difference in 
slit PA. for both datasets is relatively small (23°) and the spatial 
profiles for the WHT spectrum (not shown here) are very similar 
to those of the NOT (Fig. [2]i a likely consequence of the round 
symmetry of the object. On the other hand, in the WHT data the 
slit passes through the star partly overlapping the nebula on its 
eastern side (see Fig. [TJ, and special care was taken to extract 
this star from the nebular spectrum. 

A one-dimensional spectrum of the object (both from the 
WHT and the NOT data) was extracted summing the nebular 



spectra along the spatial axis. The WHT spectrum is shown in 
Fig.El 

Line fluxes were measured fitting Gaussian profiles using 
the task splot of IRAF. The errors in the line fluxes were cal- 
culated considering both the statistical errors and errors due to 
the flux calibration. The latter were obtained for each line from 
the error in the calibration data at the wavelength under ques- 
tion and the local RMS of the fitted calibration curve, and were 
added quadratically to the statistical errors. Resulting line fluxes 
are listed in Tables Q] (WHT) and [2] (NOT); fluxes in common 
in both datasets were found to agree within 3-cr, and in most 
cases within 1-cr. The WHT line fluxes were used in the physico- 
chemical analysis that follows. 

3.1. Reddening 



The reddening law of iFitzpatrickl (12004 for R v = 3.1, and the 
B aimer line fluxes from Tables \T\ and [2] were used to derive the 
extinction coefficient Cp from both the WHT and NOT spectra. 
For the WHT data, we derive c Hp = 0.7 + 0.1 and 0.9 + 0.1 from 
the Hy/H/3 and H5/H/3 line ratios, respectively. For the NOT, very 
similar values are obtained, namely cnp = 0.8+0.1, 0.7+0.2, and 
0.6+0.2 from the Ha/H/3, Hy/HyS, and H6/Hfi ratios, respectively. 
The weighted mean of all measurements is cnp = 0.8 + 0.1, i.e. 
A y = 1.7 ±0.2. 

We have adopted this value to deredden the NOT spectrum 
and, separately, the blue and red WHT spectra of the nebula. 
Given the lack of a precise match of the blue and red parts in the 
WHT data (Sect. 2.4) a final correction to these data was applied 
after dereddening by rescaling the whole red spectrum to m atch 
the theoretical value /(Ha)//(H/3) = 2.86 dOsterbrockll 19741) . 

The extinction derived in this paper is within the errors, al - 
beit slightly larger, than the value from Mampa so et al.l (120051) . 
Ay = 1 .5 +0.2 mag (derived from a noiser HafH/3 ratio and using 
a different extinction law). 

3.2. Density and temperature 

The IRAF package temden was used to derive the electronic den- 
sity and temperature. The density was derived by iterating the 
density calculated from the [Sn]6716 A/6731 A line ratio and 
the temperature derived from the [Nn](6548 A+6583 A)/5755 
A line ratio. This leads to N e =190+120 cirT 3 , when the errors 
in the dereddened sulphur line fluxes are considered. The density 
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Fig. 3. The WHT+ISIS spectrum of IACPN. 



error is large as the [S n]6716 A/6731 A line ratio is at the upper 
limit of its usability for the N e determination. 

Adopting the [S n] density, the electronic temperature 
was derived from [Nn](6548 A+6583 A)/5755 A, [On](3726 
A+3729 A)/(7320 A+7330 A), and [Om](4959 A+5007 
A)/4363 A line ratios giving, respectively, 11700+800 K, 
12000+800 K, and 11700+300 K. 



3.3. Chemical abundances 

The ionic abundances of He + and He ++ a nd correspondin g errors 
were calculated using the equations of iBeniamin et al.1 (1999) 
which include the effects of collisional excitation on Hei lines. 
The ionic abundances for the rest of the elements were cal- 
culated using the nebu lar analysis package in IRAF/STSDAS 
(IShaw & Dufouri 11994 . For the ions with the lowest ionisa- 
tion potentials, + ([On] 3727 A, 7319 A, 7330 A), N + ([Nn] 
5755 A, 6548 A, 6583 A), and S + ([S n] 6716 A+ 6731 A) the 
r c [Nn] was adopted, while for the higher ionisation potential 
ions, He + (He i 4471 A, 5876 A, 6678 A), He ++ (Hen 4686 A), 
++ ([Om] 4363 A, 4959 A, 5007 A), Ne ++ ([Nem] 3869 A), 
S ++ ([S m] 6312 A, 9069 A), and Ar ++ ([Arm] 7135 A, 7751 A), 
the r e [Om] was used. The [Sn] density was adopted for all the 
ions. The errors were calculated using the Monte Carlo method, 
assigning for each variable (flux, temperature, and density) 100 
random values within a normal distribution defined by the cor- 
responding quantity and its 1-cr uncertainty and calculating the 



resulting RMS in the ionic abundance. The ionic abundances de- 
rived from different lines of the same ion agree within 1 .5cr in all 
cases. The total elemental abundances were then determined by 
applying ionis ation correction factors (ICF s) following the pre- 
scriptions by Kingsburg h & Barlowl (Qj94), and the errors were 
propagated to include the ionic abundance errors derived above 
and the errors in the ICF's. The ionic and total abundances with 
their corresponding errors are given in Table 



4. Distance 

The distance to IACPN was estimated using two different meth- 
ods. 

Firstly, we used the extinction-distance t echnique (A , - D \ 
c.f. iGafhier et all Il986l) as implemented by ISale et al.1 (2009) 
for IPHAS data. The A v - D method is currently the only one 
which is capable of determining accurate (~ 30%) distances for 
a large number of PNe, being independent of statistical assump- 
tions about the physics o f the nebulae and/or of their central stars 
(Gia mmanco et alJ uOl 1). GAIA will hopefully improve this sit- 
uation in the near future by providing hundreds of trigonometric 
PN distances. The A,, - D curve towards the sightline of IACPN 
(/, b = 178.13, -4.04) was calculated using photometry of 3629 
stars in a 10' x 10' box centred on the position of the PN and 
is shown in Fig. [4] Most of the interstellar extinction builds up 
within the first kiloparsec, beyond which it remains roughly con- 
stant at Ay ~ 2. We note that this behaviour is consistent with 
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Table 1. WHT+ISIS observed line fluxes F (Obs.) normalised 
to F(H/3) = 100 in the blue section of the spectrum (separated 
by a horizontal line) and to F(Har) = 100 in the red section. See 
text for details about the normalisation of the dereddened fluxes, 
F (deredd). The wavelengths are in units of A and the total er- 
rors in the flux measurements and the propagated errors in the 
dereddened flux are given within brackets in percentage. 
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the extinction map of ISchlegel etal] d!998l) which gives an inte- 
grated Galactic extinction of Ay =1.9 towards this sightline, i.e. 
close to the asymptotic value shown by Fig. |4] The likely rea- 
son for such an abrupt asymptoting is that, at the distance of the 
Perseus arm (the only major dusty structure encountered beyond 
1 kpc along this line of sight), a Galactic latitude of b — -4.04° 
corresponds to a distance of ~ 140 pc below the midplane, and 
little cumulative dust obscuration is expected farther out. 

According to Fig. [4] IACPN could be located at any distance 
from <1 kpc (the formal best fit to the uprising section of the 
curve) to >8 kpc (the limit where there are adequate stellar data). 
The error bars of the stellar and the PN data overlap at 2cr and the 
plateau of interstellar extinction along this line of sight precludes 
a thrustworthy estimation of the extinction distance to this PN. 

Secondly, the relationshi p between Hg-surface brig htness 
and nebular radius (S - r) of iFrew & Parker! (l2006[ [2010) was 
used. IACPN is almost circular, and its dimensions, measured 
directly from the Ha IPHAS image at the 10 percent of peak- 



Table 2. As in Table 1 for the NOT spectrum normalised to 
F(H/?)=100. 
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5754 


.8 


4.4 (21) 


3.0 (22) 


Hei 5875.64 


5874 


.8 


22.1 (5) 


14.5 (7) 


[Oi] 6300.30 


6299 


.8 


41.7 (5) 


23.7 (9) 


[S in] 6312.10 


6310 


.6 


2.5 (57) 


1.4 (58) 


[O i] 6363.78 


6362 


.7 


14.5 (7) 


8.1(11) 


[Nn] 6548.03 


6546 


.9 


96.5 (4) 


51.0 (9) 


Ha 6562.82 


6561 


.8 


553.9 (4) 


291.3 (9) 


[Nn] 6583.41 


6582 


.5 


311.4 (4) 


162.7 (9) 


Hei 6678.15 


6677 


.3 


8.4(18) 


4.3 (20) 


[Sn] 6716.47 


6715. 


.2 


37.7 (5) 


19.0(10) 


[S ii] 6730.85 


6729 


.5 


29.4 (5) 


14.7 (11) 



Table 3. Ionic and total abundances for IACPN. The percentage 
errors are given within brackets. 



Ion/Element 


Abundance 


12+log(X/H) 


Hei 


0.117(10) 




Hen 


0.009 (7) 




He/H 


0.125 (9) 


11.10±0.04 


[On] 


7.69E-5 (4) 




[Om] 


1.42E-4 (2) 




icfiO) 


1.05 (8) 




O/H 


2.30E-4 (7) 


8.36±0.03 


[Nn] 


2.07E-05 (0.2) 




icm 


2.99 (8) 




N/H 


6.20E-05 (8) 


7.79±0.03 


[Ar m] 


5.72E-07 (1) 




ieKM) 


1.87 (22) 




Ar/H 


1.07E-06 (22) 


6.03±0.10 


[Sn] 


6.14e-07 (18) 




[Sn.] 


1.50E-06 (3) 




icf(S) 


1.12(1) 




S/H 


2.37E-06 (6) 


6.37±0.02 


[Nem] 


5.26E-05 (12) 




icf(m 


1.62 (7) 




Ne/H 


8.50E-05 (14) 


7.93±0.06 



bright ness contour, following the procedure of iTvlenda et alj 
(2003), are 10.5 x 9.7 arcsec. 

The total Ha flux was calculated from the wide-slit H/3 flux 
(Section 2.2) and the Balmer decrement of 5.5 + 0.2, resulting in 
log F(Ha) = -12.73 + 0.05. The Ha flux was also independently 
estimated from the IPHAS calibrated Ha image after subtracting 
the contribution of the two adjacent [N n] lines as measured form 
the spectrum, yielding log F(Ha) = -12.96 + 0.22, in fair agree- 
ment with the spectroscopically determined value. We will adopt 
the better determined spectroscopic value in the following. 

The resulting Ha surface brightness is log S(Ha) = -4.00 + 
0.06 erg cm -2 s~' sr~'. After dereddening this value using Ay = 
1.7 + 0.2 mag and the extinction law of Fitzpatrick (2004) with 
Rv = 3.1, we find the physical radi us of the PN to be 0. 31 pc, 
using the mean Ha S — r relation of lFrew & Parker! d2006). The 
mean-trend equation (Frew 2008, Frew et al. 201 1, in prepara- 
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Fig. 4. IPHAS extinction - distance curve (circles) for field stars 
towards the sightline of IACPN. The horizontal line shows the 
measured extinction Ay of IACPN flanked by two dashed lines 
representing ±l<x error limits. 

tion) was used, which is applicable to round, optically-thick PN, 
viz.: 

logS(Ha) = -3.65(+0.08)logfl - 5.34(+0.07). (1) 

The resulting distance is 12.8 + 3.7 kpc. The quoted uncer- 
tainty in the di stance inclu des the formal dispersion of the re- 
lation (+ 28%: iFrewl 120081) which is the dominant error term, 
added in quadrature with the estimated uncertainties of the in- 
tegrated Ha flux, mean diameter, and reddening of the PN. The 
corresponding galacto-centric distance i s Pgr _ 20.8 + 3 .8 kpc 
if a distance to the Sun of 8.0 + 0.6 kpc dGhez et al.ll2008l) is as- 
sumed. At 12.8 kpc from the Sun, IACPN is located about 1 kpc 
off the plane of the Galactic disk. That is a large height indeed, 
but we note that th e well-known flaring of the outer disk (see e.g. 
iGvuk et alii 19991) allows a high probability of thin disk objects 
being located eve n at this offset from the mid-plane. Using the 
Gvuk et alj (1 19991) prescription, the thin disk scale height at Dec 
= 20.8 kpc is in excess of 1 kpc (even i f the solar neigh bourhood 
scale height is reduced in line with the Juric et al ] (120081) result). 

We will adopt for IACPN the S - r distance of 12.8 kpc. 
At Dcc=20.8 kpc, it would be one of the most distant planetary 
nebula from the Galactic Centre for which a distance has been 
determined (c.f. lAcker et al.lll994t IStanghellini et al.ll2008l) and 
the fart hest PN from the Galactic Centre with m easured abun- 
dances dMaciel & Costdl2009fclHenrv et al.ll2010l) . 

5. Discussion 

IACPN is a low density planetary nebula with an emission line 
spectrum typical of a moderate e xcitation Type II PN (Class 3 of 
Dopita & Meatheringham| [r990l) . This, together with its strong 
stratification, apparent from Hen to Oi (Fig.|2]l, suggests that it 
is composed of a moderately massive nebula hosting a relatively 
hot star. The ionised hydrogen mass of the nebula, derived from 
its H/3 flux, electron temperature, radius, and distance, is M= 0.4 
M . The nebular expansion velocity (17 km s~') and the physical 
radius (0.31 pc at the distance of 12.8 kpc) lead to a kinematic 
age of 17800 yr. Such a large kinematical age would point to a 
low mass centr al star if a lumi nosity larger than ~ 300 - 500 L 
is required (cf. lBlockerlll995l) . To crud ely estimate the ce ntral 
star mass, we follow the prescriptions of M aciel et al.l (1201 Ol) for 



a N/O abundance ratio as given in Table [3] finding m cs = 0.6 
Mo and an age for the progenitor star of ~4 Gyr (if Case A from 
iMaciel et alj|2010l is assumed). However, given the several un- 
certain parameters involved (central star luminosity and temper- 
ature, kinematical vs. real age), no precise estimate of the stellar 
mass can be derived. 

We have ran several simple Cloudy dFerland et al.l Il998l) 
photoionisation models usin g spherical or shell geometry and 
the model atmospheres of Rauch (2003). We assumed a nebula 
with an outer radius of 0.31 pc, hydrogen density of 200 cm -3 , 
and chemical abundances 0.3 dex below those labelled in Cloudy 
as "planetary nebula". This implies that the input models have 
12+log(0/H)=8.34 and He/H=0.1, close to our measured val- 
ues. On the other hand, the model abundances of N, Ar, and S are 
0.1-0.3 dex above our derived values, whereas the model abun- 
dance of Ne is 0.2 dex below the value listed in Table|3] Since the 
abundances of these elements are more uncertain, and since our 
comparison is based on the He and O ions, this should be a good 
approximation for our purposes. We considered the observed in- 
tensities of the He i, He n, [O i], [O n], and [O m] lines, as given in 
Table [1] and found that high stellar temperatures (T e ff Z 10 5 K) 
are needed to account for the observed intensity of the He n 4686 
A line, in agreement with the results presented above. Models 
with those effective temperatures also reproduce the intensities 
of all the aforementioned lines. The radiation-bounded case re- 
quires a luminosity for the ionising source L CJ ~ 100 L Q , being 
therefore consistent with a low mass star on the cooling track. As 
a matter of fact, similar models for radiation-bounded nebulae as 
small as 0.01 pc (corresponding to extinction distances smaller 
than 1 kpc; see Section 4) would imply central star luminosities 
around 0.001 L Q , i.e., much lower than the expected luminosity 
range of 10-1000 L for stars on the cooling track. This renders 
implausible models using very short distances for this PN. 

We also investigated the [O in] luminosity of IACPN to com- 
pare it to the bright-end cutoff of the PN luminosity function 
(PNLF: ICiardullol l20Toh . We used the long-slit integrated H/3 
flux and the [Om]/Hy6 ratio of 7.3 to derive an integrated [O in] 
flux, F(5007) = 2.43 x 10~ 13 erg cm ' 2 s" 1 . Correcting for ex- 
tinction, and following Jacoby ( 1989), we derive a PN absolute 
magnitude, M5007 = +0.3*? ! at the adopted distance. The PN is 
thus ~4.8 magnitudes down the PNLF from the bright cut-off 
magnitude, M* = -4.48. This is consistent with the evolved na- 
ture of the PN, and the position of the (unobserved) central star 
on the WD cooling track. 

In summary, all physical data at hand on IACPN are con- 
sistent with a distant, aged, radiation-bounded Type II planetary 
nebula excited by a low mass central star with relatively high 
temperature and having already reached its cooling track. 

With respect to the chemistry, the O/H abundance of 
8.36+0.03 (Table |3]l is interestingly high for the large galacto- 
centric distance, 20.8 kpc, and makes IACPN a key target for 
studies of the abundance gradient in the Milky Way. It is gen- 
erally agreed that there is some negative PNe metallicity gradi- 
ent but the consensus seems to finish there. Values for the best 
measured element, oxygen, range from Alog(Q/ H)/D(7c = -0.01 
dPottasch & Bernard-Salasl 120061) to -0.08 ( Stanghellini et al.l 
2006) with most recent values still disagreeing by a factor of two, 
well above the uncertainties qu oted by the authors: cf. -0.023 + 
0.00 6 and -0.058 + 0.00 6 from lStanghelhni & Havwoodl (l2010l) 
and iHenrv et"al] (|2010), respectively. Similar confusion ham- 
pers the possible existence of a variable slope, i.e. if the 
gradient steepen s (Stanghel lini & Havwoodl 1201 Oh or flattens 
(Ma ciel & Costal 1201 Ol) outwards along the Galactic equato- 
rial plane, and also on its secular temporal behaviour, either 
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null (lHenrvetal.1 12010) or ste eping (IStanghellini & Havwoodl 
1201 $i or flattening with time dMaciel & Costal l2010l). Perhap s 
the most sensible conclusion comes from iHenrv et al.l (|2010): 
"Essentially we have reached a confusion limit regarding the 
abundance gradient as derived using PNe.". One of the major 
sources of this confusion is the unsettled PNe distance deter- 
mination problem allied with the extreme paucity of nebulae at 
large galac tocentric dist a nces. I n fact, in the most recent works, 
those by IHenrv et all (1201 Oh and IStanghellini & Havwoodl 
(l2010l) . only two PNe with Dec > 17 kpc are included in each 
sample (K3-64, K3-70, and Ml-9, K3-66, respectively), all four 
PNe having mediocrily determined O abundances (quoted errors 
from 0.2 to 0.5 dex) and distances. However, these distant nebu- 
lae are obviously critical for the gradient determination. 

In this context, IACPN provides a valuable data point located 
at the outer reaches of the Galaxy, and its well-constrained O 
abundance do support models with a flatter gradient toward the 
external disk. An intriguing possibility is that IACPN and per- 
haps also the other four PNe with Dec >17 kpc might even have 
an extragalactic origin, if indee d the Monoceros Ring proposed 
to be located at these radii (eg. IConn et al]|2008l) is attributable 
to a satellite accret ion event. A consensus on this has no t yet 
been reached (eg. Hammerslev & Lop ez-Corredoirall20 1 Ob . But 
either way it is clear that the properties of the outer Galactic disc 
are still in need of much better observational definition. 

6. Conclusions 

We are exploring the new IPHAS planetary nebula candidates 
dViironen et al.ll20"09al) towards the Galactic Anticentre in order 
to better sample the Galactic abundance gradient at large galac- 
tocentric distances. In this paper, the discovery of a new PN lo- 
cated at only 4 degrees from the Galactic Anticentre, IPHASX 
J052531. 19+281945.1, is presented, its physical and chemical 
properties are studied and its distance determined. The object 
turned out to be a low-density, distant, and aged PN. It is the 
Galactic PN with the largest galactocentric distance (20.8+3.8 
kpc) where the chemical abundances have been measured, and 
this, combined with its relatively high oxygen abundance, 12 
+ log(0/H) = 8.36+0.03, points to a flattening of the Galactic 
abundance gradient towards the outer Gal axy rather than to a lin- 
early decreasing oxygen abundance (e.g. lMaciel & Costal l2009, 
and references therein). 

To properly study the behaviour of the abundance gradient 
at the outer Galaxy, reasonably accurate abundance measure- 
ments on a sample of objects at large galactocentric distances are 
needed. We have recently discovered in the IPHAS images about 
180 new faint PN candid ates located in a 30° x 10° region around 
the Anticentre dViironen et al.l l"2009a). Further abundance mea- 
surements of confirmed candidates in this region would help to 
build a robust chemical gradient determination and reach firmer 
conclusions concerning its flattening. 
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